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ABSTRACT 


The  structure  of  the  heliospheric  current 
sheet  on  a  spherical  source  surface  of  radius  2.35 
R0  has  been  computed  using  a  potential  field  model 
during  the  first  year  and  a  half  after  the  last 
sunspot  minimum.  The  solar  polar  magnetic  field 
that  is  not  fully  observed  in  conventional  magne¬ 
tograph  scans  was  included  in  the  computation. 
The  computed  heliospheric  current  sheet  had  a 
guasi-stationary  structure  consisting  of  two 
northward  and  two  southward  maxima  in  latitude  per 
solar  rotation.  The  extent  in  latitude  slowly 
increased  from  about  15  degrees  near  the  start  of 
the  interval  to  about  45  degrees  near  the  end  of 
the  interval.  The  magnetic  field  polarity  (away 
from  the  Sun  or  toward  the  Sun)  at  the  sub¬ 
terrestrial  latitude  on  the  source  surface  agreed 
with  the  interplanetary  magnetic  field  polarity 
observed  or  inferred  at  Earth  on  82%  of  the  days. 
The  interplanetary  field  structure  observed  at 
Barth  at  this  time  is  finely  tuned  to  the  struc¬ 
ture  of  low-latitude  fields  on  the  source  surface. 
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i*  Introduction 

Almost  daily  observations  of  the  large-scale  photos- 
pheric  magnetic  field  structure  were  started  at  the  Stanford 
Solar  Observatory  in  May  1976  and  have  continued  to  the 
present  time.  We  compute  the  large-scale  structure  of  the 
magnetic  field  in  the  heliosphere  using  Zeeman  observations 
of  the  1 ine-of-sight  component  of  the  photospheric  magnetic 
field  together  with  a  potential  field  model.  It  is  also 
possible  to  infer  the  structure  of  the  heliospheric  current 
sheet  from  the  maximum  brightness  contours  in  the  K  coronam- 
eter  observations  at  Mauna  Loa  Observatory  (Burlaga  et  al. 
1981  and  references  therein) . 

During  the  time  interval  considered  here  there  was  an 
electric  current  sheet  that  was  warped  northward  and  south¬ 
ward  of  the  plane  of  the  solar  equator  (Schulz,  1973). 
North  of  the  current  sheet  the  interplanetary  magnetic  field 
(IMP)  was  directed  away  from  the  Sun  and  south  of  the 
current  sheet  the  IMP  was  directed  toward  the  Sun 

The  minimum  between  sunspot  cycles  20  and  21  occurred 
in  June  1976.  During  the  18  Carrington  Solar  Rotations 
beginning  in  May  1976  the  computed  current  sheet  was  quasi¬ 
stationary,  having  in  each  solar  rotation  two  northward 
extensions  and  two  southward  extensions.  This  usually 
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produced  the  characteristic  four  sector  structure  in  the 
interplanetary  magnetic  field  observed  at  Earth  (Svalgaard 
and  Wilcox,  1975).  Occasionally  during  a  rotation  one  or 
even  both  of  the  northward  extensions  of  the  current  sheet 
"missed”  the  Earth  resulting  in  a  two  sector  or  even  a 
"zero”  sector  structure  being  observed  at  Earth.  Around 
sunspot  minimum  the  maximum  extent  in  latitude  of  the  com¬ 
puted  current  sheet  was  about  15  degrees,  while  by  the  end 
of  the  18  solar  rotations  discussed  here  the  maximum  lati¬ 
tude  had  increased  to  about  45  degrees.  -Just  after  the  time 
interval  discussed  here  the  maximum  latitude  of  the  current 
sheet  increased  further  and  the  quasi-stationary  structure 
of  the  current  sheet  began  to  change,  so  that  September  1977 
seems  a  natural  point  to  end  the  present  investigation.  The 
structure  of  the  computed  heliospheric  current  sh?et  in 
later  portions  of  sun3pot  cycle  21  will  be  discussed  in 
future  papers. 

We  also  investigate  here  the  effect  of  varying  the 
source  surface  radius,  the  strength  of  the  polar  field 
correction  and  the  latitude  on  the  source  surface  used  to 
predict  the  IMF  polarity  seen  at  Earth. 

2.  Computation  of  the  Source  Surface 

Schatten  et  al.  (1969)  and  Altschuler  and  Newkirk 


(1969)  introduced  the  concept  of  a  potential  field  model 
with  a  spherical  source  surface  surrounding  and  concentric 
with  the  Sun  (see  also  Levine  and  Altschuler,  1974;  Poletto 
et  al.,  1975;  Altschuler  et  al.,  1976;  Adams  and  Pneuman, 
1976;  Svalgaard  and  Wilcox,  1978  and  Riesebieter  and  Neu- 
bauer,  1979).  Outside  the  source  surface  it  is  assumed  that 
the  radial  flow  of  the  solar  wind  carries  the  magnetic  field 
outward  into  the  heliosphere.  Between  the  photosphere  and 
the  source  surface  it  is  assumed  that  the  magnetic  field  can 
be  described  in  terms  of  a  potential  which  satisfies 
Laplace's  equation.  For  the  work  described  here  the  inner 
boundary  condition  at  the  photosphere  is  the  1 ine-of-sight 
magnetic  field  observed  at  the  Stanford  Solar  Observatory. 
The  outer  boundary  condition  is  that  the  field  is  normal  to 
the  source  surface,  consistent  with  the  assumption  that  it 
is  then  carried  outward  by  the  solar  wind.  The  assumption 
in  the  3ource  surface  calculation  that  there  are  no  currents 
would  not  be  very  good  for  the  3trong  localized  fields  of  an 
active  region,  but  for  the  large-scale,  quasi-stationary 
fields  that  dominate  the  present  analysis  the  source  surface 
gives  a  reasonably  good  prediction  of  the  polarity  of  the 
interplanetary  magnetic  field  observed  at  Earth. 

A  non-spher ical  source  surface  computation  (Schulz  et 
al. ,  1978;  Levine  et  al., 


1982)  should  give  an  improved 
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prediction  of  the  coronal  structure  and  of  the  IMP  observed 
at  Barth,  but  as  we  shall  see  the  spherical  source  surface 
already  does  quite  well  at  predicting  the  IMP  polarity.  The 
amount  of  improvement  to  be  obtained  from  a  non-spherical 
source  surface  using  our  observations  will  be  investigated 
in  a  later  paper.  The  magnitude  of  the  interplanetary  mag¬ 
netic  field  may  be  better  corouted  with  non-spherical  source 
surface . 

In  most  previous  work  the  magnetic  field  on  the  source 
surface  has  been  computed  only  once  for  each  Carrington 
Rotation,  i.e.  ii  steps  of  360  degrees  in  longitude.  This 
forces  the  beginning  (360  degrees)  and  the  end  (0  degrees) 
of  the  rotation  to  have  the  same  structure,  even  though  they 
are  separated  in  tine  by  27  days.  To  avoid  this  difficulty 
we  have  computed  the  field  on  the  source  surface  in  steps  of 
10  degrees  in  tha  starting  longitude,  and  retained  only  the 
central  interval  of  width  30  degrees  in  longitude  from  each 
such  computation.  As  the  last  step  a  1:2:1  averaging  of  the 
three  calculations  for  each  longitude  strip  is  applied  to 
slightly  smooth  the  structure.  In  the  (rare)  case  of  data 
gaps  we  interpolate  between  the  previous  and  the  subsequent 
rotation . 

Stenflo  (1971),  Howard  (1977),  Svalgaard  et  al.  (1978), 
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Pneuraan  et  al.  (1978),  and  Burlaga  et  al.  (1981)  have 
pointed  out  that  conventional  solar  magnetograph  observa¬ 
tions  do  not  adequately  represent  the  solar  polar  magnetic 
£ield  strength.  Wilcox  et  al.  (1980)  computed  the  helios¬ 
pheric  current  sheet  configuration  early  in  1976  using  solar 
magnetograph  observations  from  Mt.  Wilson  Observatory  that 
were  not  corrected  for  the  solar  polar  magnetic  field  not 
observed  in  daily  solar  magnetograph  observations.  As  a 
result  the  computed  extent  in  latitude  of  the  heliospheric 
current  sheet  was  probably  too  large,  as  was  pointed  out  by 
Burlaga  et  al .  (1981). 

In  the  present  computation  of  the  heliospheric  current 
sheet  we  use  the  magnetic  field  observed  at  a  resolution  of 
three  arc  minutes  in  daily  scans  with  the  solar  magnetograph 
at  the  Stanford  Solar  Observatory,  plus  the  solar  polar 
field  strength  determined  by  Svalgaard  et  al.  (1978)  for  the 
same  solar  rotations  analyzed  in  the  present  paper.  In  the 
interval  analyzed  by  Svalgaard  et  al.  (1978)  the  magnitude 
of  the  solar  polar  field  did  not  change  appreciably.  We 
note  that  near  the  minimum  of  the  sunspot  cycle  the  solar 
polar  fields  will  have  the  maximum  influence.  As  the  sun¬ 
spot  cycle  progresses  after  minimum  the  strength  of  the 
solar  polar  field  decreases  while  the  strength  of  the  low 
latitude  fields  increases.  Near  sunspot  maximum,  when  the 


The  radial  magnetic  field  computed  on  a  spherical 
source  surface  at  2.35  for  Carrington  Rotation  1648 
beginning  7  November  1976  is  shown  in  Figure  1.  The  current 
sheet  is  represented  by  the  zero  contour  shown  as  a  thick 
solid  line  near  the  equator.  The  solid  contours  above  it 
represent  field  directed  away  from  the  Sun  with  relative 
magnitudes  1,  5,  and  10,  while  the  dashed  contours  represent 
field  directed  toward  the  sun  with  the  same  contour  levels. 
The  predominance  of  away  polarity  magnetic  field  in  most  of 
the  northern  region  of  the  heliosphere  and  of  toward  field 
in  most  of  the  southern  heliosphere  is  apparent  in  Figure  1. 

The  +  (away  from  the  Sun)  and  -  (toward  the  Sun)  sym¬ 
bols  in  Figure  1  represent  daily  polarities  of  the  inter¬ 
planetary  magnetic  field  at  Earth  as  observed  by  spacecraft 
(King,  1979)  or,  when  spacecraft  observations  were  not 
available,  inferred  from  polar  geomagnetic  observations 
(Svalgaard,  1973).  The  INF  polarities  at  Earth  plotted  in 
Figure  1  have  been  displaced  by  five  days  corresponding  to 
the  average  transit  time  of  solar  wind  from  Sun  to  Earth 


3 


near  the  times  when  the  large-scale  magnetic  polarity 
changes  (sector  boundaries).  The  average  magnitude  of  this 
transit  time  during  the  solar  rotations  studied  here  has 
been  determined  by  a  cross  correlation  analysis  to  be 
described  later.  We  note  that  near  the  sector  boundaries 
the  velocity  of  the  solar  wind  is  almost  always  near  a 
minimum  (Wilcox  and  Ness,  1965),  so  that  this  transit  time 
is  longer  than  the  average  solar  wind  transit  time. 

Figure  2,  in  the  same  format  as  Figure  1,  shows  the 
field  computed  at  the  source  surface  for  Carrington  Rotation 
1656  beginning  13  June  1977.  The  extent  in  latitude  of  the 
computed  current  sheet  had  increased  to  about  40  degrees, 
but  the  same  property  of  two  northward  excursions  and  two 
southward  excursions  in  the  current  sheet  (a  four  sector 
structure)  was  still  evident. 

Figures  3a  and  3b  show  the  computed  currant  sheets  and 
IMF  polarities  observed  at  Earth  during  the  18  solar  rota¬ 
tions  considered  in  the  present  work.  In  every  rotation 
except  Number  1644  there  were  two  northern  and  two  southern 
extensions  of  the  currant  sheet,  corresponding  to  a  basic 
four  sector  structure.  In  Rotation  1645  the  computed 


current  sheet  was  everywhere  southward  of  the  heliographic 
latitude  of  the  Earth,  and  the  IMF  polarity  observed  at 


Earth  was  almost  entirely  away  from  the  Sun.  This  presum¬ 
ably  is  an  example  of  the  situation  discussed  by  Wilcox 
(1972)  in  which  near  the  last  five  (now  six)  sunspot  minima 
the  observed  or  inferred  IMF  polarity  has  been  largely  away 
from  the  Sun  during  a  few  consecutive  rotations.  If  the 
current  sheet  "misses"  the  Earth  near  the  time  of  a  sunspot 
minimum  the  resulting  predominant  polarity  of  the  IMF  could 
oe  either  away  from  or  toward  the  Sun  according  to  the  con¬ 
siderations  discussed  in  this  paper.  A  predominance  in  away 
polarity  in  the  observed  photospheric  field  also  discussed 
by  Wilcox  (1972)  would  not  necessarily  be  directly  related 
to  the  situation  shown  here  in  Rotation  1645. 

Hundhausen  (1977)  noted  that  a  "monopolar"  sector 
structure  as  seen  in  Rotation  1645  of  Figure  3a  might  appear 
at  the  beginning  of  a  new  solar  cycle.  However,  the  sugges¬ 
tions  that  at  this  time  "The  prominent  recurrent  sectors, 
streams  and  geomagnetic  activity  sequences  should  end 
abruptly"  and  the  "Recurrence  with  the  27-day  solar  rotation 
period  should  become  rare"  are  not  consistent  with  the  com¬ 
puted  current  sheets  in  Figures  3a  and  3b. 

In  Rotation  1658  the  computed  current  sheet  had  a  clear 
"four  sector"  structure,  but  was  sufficiently  far  south  of* 
the  heliographic  latitude  of  the  Earth  that  only  a  two  sec¬ 
tor  structure  was  observed  here.  This  appears  to  be  the 
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same  geometry  but  the  opposite  sense  from  the  situation  in 
early  1976  described  by  Scherrer  et  al.  (1977). 

From  the  start  of  Figure  3a  near  the  minimum  of  the 
eleven  year  sunspot  cycle  to  the  end  of  Figure  3b,  1.5  years 
later,  the  maximum  extent  in  latitude  of  the  computed 
current  sheet  increased  from  about  15  degrees  to  about  45 
degrees.  This  increase  is  qualitatively  similar  to  but 
larger  than  the  average  variation  computed  by  Svalgaard  and 
Wilcox  (1976)  through  the  previous  four  sunspot  cycles. 

Burlaga  et  al .  (1981)  noted  that  for  Carrington  Rota¬ 
tions  1639  and  1640,  just  before  the  start  of  the  interval 
shown  in  Figure  3a,  a  solar  dipole  magnetic  axis  tilted 
about  20°  to  15°  with  respect  to  the  solar  rotation  axis 
cannot  explain  the  sector  pattern  observed  by  Helios.  The 
sector  patterns  shown  in  Figure  3a  and  3b  during  1.5  years 
after  the  rotations  discussed  by  Burlaga  et  al.  (1981)  also 
cannot  be  explained  with  a  tilted  dipole,  as  was  proposed  by 
Smith  and  Tsurutani  (1979),  Villante  et  al .  (1979),  Smith 

and  Wolfe  (1979),  Zhao  and  Hundhausen  (1981),  and  Hakamada  and  Akasofu 
(1981). 

On  most  of  the  rotations  during  1976  shown  in  Figure  3a 
the  current  sheet  extended  more  into  the  southern  helio¬ 
sphere  (the  case  of  Rotation  1644  is  discussed  below),  con¬ 
sistent  with  the  results  of  Wilcox  et  al  .(1980),  Burlaga  et  al . 
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(1981)  and  Villante  et  al .  (1982).  The  conjecture  o£  Vil- 
lante  et  al .  (1982)  that  the  current  sheet  during  the  first 
half  of  1977  was  confined  in  a  narrower  latitude  region  is 
not  consistent  with  the  current  sheets  shown  in  Figure  3b. 

In  Figures  3a  and  3b  intervals  of  significant  disagree¬ 
ment  between  the  IMF  polarity  predicted  by  the  computed 
current  sheet  and  that  actually  observed  are  indicated  by  a 
bar  attached  to  the  current  sheet.  We  note  that  for  the 
most  part  the  daily  polarity  of  the  IMF  observed  at  Earth  is 

quite  well  predicted  by  the  computed  current  sheet;  in  fact 
* 

there  is  agreement  on  32%  of  the  days. 

A  conspicuous  disagreement  is  associated  with  the  rapid 
change  in  the  computed  current  sheet  from  one  rotation  to 
the  next  at  Rotation  1644.  This  change  in  the  computed 
current  sheet  was  caused  by  the  appearance  of  a  particularly 
large  bipolar  magnetic  region  in  the  photosphere.  Figures 
4a,  4b,  and  4c  show  synoptic  charts  of  the  observed  photos- 
pheric  magnetic  field  for  Rotations  1643,  16*4  and  1645.  A 
large  bipolar  magnetic  region  appeared  in  Rotation  1644  at 
longitude  120  degrees  with  predominantly  toward  polarity 
field.  The  corresponding  IMF  polarity  observed  at  Earth  was 
away  on  several  days  during  which  the  computed  current  sheet 
would  lead  to  a  prediction  of  toward.  It  seems  possible 
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that  there  may  have  been  a  region  of  toward  magnetic  field 
polarity  in  the  heliosphere  corresponding  to  this  bipolar 
magnetic  region,  but  at  a  latitude  sufficiently  far  north  so 
as  not  to  intersect  the  Earth.  A  similar  event  occurred 
near  140  degrees  longitude  in  the  southern  hemisphere  in 
Carrington  Rotation  1651. 

The  rather  rapid  change  in  the  computed  current  sheet 
near  longitude  zero  from  Rotations  1652  to  1653  was  also 
caused  by  the  appearance  of  a  large  bipolar  magnetic  region 
in  the  photosphere,  but  in  this  case  the  region  remained  in 
the  photosphere  for  several  rotations,  and  the  corresponding 
effects  on  the  computed  current  sheet  also  continued  for 
several  rotations. 

In  many  of  the  rotations  shown  in  Figures  3a  and  3b, 
the  latitude  of  the  current  sheet  at  the  end  of  the  rotation 
is  significantly  different  from  the  latitude  at  the  start  of 
the  rotation.  This  illustrates  the  advantage  gained  from 
computing  the  field  structure  on  the  source  surface  at  steps 
of  10  degrees  in  the  starting  longitude,  since  if  only  one 
computation  were  made  for  each  rotation  the  latitude  of  the 
current  sheet  at  the  start  and  the  end  of  the  rotation  would 
be  forced  to  be  the  same. 
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4.  Influence  on  tbs  Solar  Polar  Field  Strength  and  the 
Radius  of  the  Source  Surface 

The  source  surface  current  sheet  was  computed  in  the 
above  discussions  using  the  solar  polar  field  strength  of 

g 

the  form  11.5  cos  8  gauss  derived  by  Svalgaard  et  al. 
(1978),  where  8  is  the  colatitude.  We  will  now  investigate 
the  effect  of  changing  the  magnitude  of  the  derived  solar 
polar  field  and  of  changing  the  radius  of  the  source  sur¬ 
face.  For  this  purpose  we  will  compute  a  cross  correlation 
between  the  IMF  polarity  predicted  from  the  computed  current 
sheet  and  that  actually  observed  or  inferred  at  Earth. 

In  order  to  determine  the  predicted  IMF  polarity  a  line 
is  drawn  on  the  source  surface  at  the  heliographic  latitude 
of  the  Earth,  i.e.  varying  from  7  degrees  north  to  7  degrees 
south  through  the  year.  This  line  is  divided  into  daily 
increments,  and  on  a  given  day  if  the  current  sheet  is 
southward  of  the  line  the  predicted  polarity  is  away,  and  if 
the  current  sheet  is  northward  of  the  line  the  predicted 
polarity  is  toward  the  Sun.  ht  least  5/8  of  a  day  must  have 
the  same  polarity  in  order  for  a  polarity  to  be  assigned. 
The  solid  curve  in  Figure  5  shows  the  maximum  cross  correla¬ 
tion  between  the  predicted  IMF  polarity  described  above  and 
the  polarity  observed  at  Earth  as  a  function  of  the  radius 


1 
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of  the  source  surface  on  which  the  current  sheet  is  corn-* 
puted,  using  the  polar  field  strength  computed  by  Svalgaard 
et  al.  (1978).  The  largest  correlation  occurs  for  a  source 
surface  of  radius  2.35  R#,  and  this  radius  has  therefore 
been  chosen  for  most  of  the  discussion  in  this  paper.  For 
comparison  Figure  5  also  shows  similar  maximum  cross  corre¬ 
lations  for  source  surfaces  computed  with  no  polar  field 
added,  and  for  5.8  gauss  and  17.3  gauss  added  polar  field. 
We  note  that  the  solar  polar  field  of  11.5  gauss  computed  by 
Svalgaard  et  al .  (1978)  does  give  the  best  agreement, 
although  the  differences  are  not  large. 

Figure  6  shows  computed  current  sheets  on  a  typical 
Carrington  rotation  (i.e.  Rotation  1656)  for  the  four  values 
of  added  solar  polar  magnetic  field.  The  current  sheet  for 
the  selected  value  of  11.5  gauss  is  shown  with  a  solid  line. 
The  current  sheet  shown  with  short  dashes  was  computed  with 
no  added  solar  polar  field,  and  it  has  the  largest  extent  in 
latitude  in  Figure  6.  The  dash  dot  line  is  a  current  sheet 
computed  with  17.3  gauss  added  solar  polar  field  (i.e.  one 
and  one  half  times  the  preferred  value) ,  and  it  has  the 
smallest  extent  in  heliographic  latitude. 

Pneuman  et  al.  (1978)  computed  the  field  on  a  source 
surface  at  2.5  R#  during  the  Skylab  period  in  1973  and  found 
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that  their  computed  neutral  lines  were  systematically  pole- 
ward  of  the  brightness  maxima  observed  at  1.8  R«  with  the 

o 

K-coronameter  at  Mauna  Loa  Hawaii.  If  the  fields  above  70 
latitude  measured  with  the  full  disk  magnetograph  at  Kitt 
Peak  National  Observatory  were  increased  to  about  30  gauss 
this  effect  was  removed.  This  is  a  much  larger  correction 
for  the  solar  polar  field  than  we  have  used.  The  reason  for 
the  difference  from  our  work  is  not  clear.  Pneuman  et  al . 
(1978)  suggested  other  possible  causes  for  their  sys¬ 
tematic  poleward  displacement  of  the  neutral  line;  to  the 
extent  that  these  operated  the  solar  polar  field  correction 
would  be  reduced.  A  difference  in  solar  magnetograph  cali¬ 
brations  between  Kitt  Peak  and  Stanford  may  have  contributed 
to  the  different  corrections,  and  the  solar  polar  field 
strength  may  have  been  different  in  1973  and  1976. 

All  the  computed  current  sheets  in  Figure  6  cross  the 
solar  equator  at  the  same  longitudes,  and  the  cross  correla¬ 
tions  shown  in  Figure  5  are  nearly  the  same  for  all  the 
values  of  added  solar  polar  magnetic  field.  Near  340 
degrees  longitude  the  maximum  latitude  of  the  current  sheet 
decreases  from  58  degrees  for  no  added  solar  polar  field  to 
37  degrees  for  17.3  gauss  added  field.  All  of  the  computed 
current  sheets  in  Figure  6  agree  almost  equally  well  with 
the  INF  polarity  observed  at  Barth. 


The  computed  current  sheet  in  Carrington  Rotation  1656 
for  several  values  of  the  radius  of  a  spherical  source  sur¬ 
face  is  shown  in  Figure  7.  As  the  radius  of  the  source  sur¬ 
face  is  increased  the  extent  in  latitude  of  the  computed 
current  sheet  decreases  since  the  relative  weight  of  the 
dipole  component  of  the  solar  magnetic  field  increases.  All 
of  the  computed  current  sheets  in  this  Rotation  agree  almost 
equally  well  with  the  observed  IMF  polarity. 

Thus  a  comparison  of  the  IMF  polarity  predicted  from  a 
computed  current  sheet  with  the  IMF  polarity  observed  at 
Earth  is  a  weak  test  of  the  extent  in  latitude  of  the  com¬ 
puted  current  sheet.  A  spacecraft  observing  at  large  helio¬ 
graphic  latitudes  would  give  the  definitive  answer  to  the 
problem  of  the  extent  in  latitude  of  the  heliographic 
current  sheet. 

5.*  Further  Comparison  of  Predicted  and  Observed  IMF  Polar- 

it* 

In  the  discussion  so  far  the  IMF  polarity  observed  at 
Earth  has  been  compared  with  the  source  surface  field  polar¬ 
ity  at  the  heliographic  latitude  of  the  Barth,  what  happens 
if  instead  we  compare  the  observed  IMF  polarity  at  Barth 
with  the  polarity  on  the  source  surface  5  degrees  north  of 
the  sub-terrestrial  latitude?  Figure  8  shows  that  the 
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maximum  cross  correlation  decreases  from  0.64  to  0.54.  We 
see  in  Figure  8  that  the  sub-terrestrial  latitude  on  the 
source  surface  has  the  most  similar  magnetic  polarity  struc¬ 
ture  to  that  observed  at  Barth,  and  that  even  a  few  degrees 
north  or  south  of  the  sub-terrestrial  latitude  the  correlation 
with  the  observed  field  is  smaller. 

For  the  adopted  conditions  of  source  surface  radius 
equal  to  2.35  Rg  and  added  solar  polar  field  of  11.5  gauss. 
Figure  9  shows  a  cross  correlation  between  the  predicted 
field  polarity  at  the  sub- terrestr ial  point  on  the  source 
surface  and  the  IMF  polarity  observed  at  Earth.  The  first 
peak  at  5.0  +  0.3  days  represents  the  transit  time  for  the 
solar  wind  plasma  to  transport  the  magnetic  field  from  Sun 
to  Earth.  The  five  day  lag  corresponds  to  a  solar  wind 
velocity  of  350  km/s.  This  represents  the  average  solar 
wind  velocity  at  sector  boundary  crossings,  which  are  usu¬ 
ally  near  minima  in  solar  wind  velocity  (Wilcox  and  Ness, 
1965).  The  relatively  slow  decline  in  amplitude  of  the 
peaks  near  32  days,  69  days  and  36  days  3hows  that  the 
large-scale  IMF  structure  is  quasi-stationary.  The  inter¬ 
mediate  peaks  are  caused  by  the  four-sector  nature  of  the 
IMF  structure  at  this  time.  The  difference  in  time  between 
the  peak  at  32  days  and  at  5  days  shows  that  the  recurrence 
time  of  the  IMF  is  close  to  27  days. 
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Our  final  comparison  of  the  structure  predicted  from 
the  source  surface  and  observed  at  Earth  is  shown  in  Pigures 
10  and  11  *  These  figures  are  now  in  a  Bartels  rotation  plot 
as  is  customary  for  geomagnetic  observations.  Comparison  of 
the  figures  shows  that  the  large-scale  structure  is  quite 
well  predicted,  and  most  of  the  disagreements  come  near  sec¬ 
tor  boundary  crossings  or  on  occasional  rotations.  A  por¬ 
tion  of  the  disagreement  near  boundary  crossings  is  caused 
by  our  use  of  a  constant  five  day  3olar  wind  transit  time 
from  Sun  to  Earth,  while  in  fact  there  are  some  variations 
among  the  actual  transit  times.  On  the  one  day  scale  used 
in  plotting  Figures  10  and  11  these  variations  in  transit 
time  would  not  be  a  large  effect. 

6.  Summary 

The  heliospheric  current  sheet  configuration  has  been 
computed  on  a  source  surface  at  2.35  R#  during  an  interval 
of  1.5  years  after  the  last  sunspot  minimum.  The  magnetic 
field  observed  on  almost-daily  scans  with  the  solar  magneto¬ 
graph  at  the  Stanford  Solar  Observatory  has  been  corrected 
for  the  solar  polar  fields  that  are  not  fully  observed. 
This  correction  significantly  reduces  the  extent  in  latitude 
of  the  computed  current  sheet. 


The  field  on  the  source  surface  at  the  sub- terrestrial 


latitude  agrees  best  with  the  interplanetary  magnetic  field 
observed  at  Earth.  A  deviation  from  this  latitude  of  more 
than  a  few  degrees  produces  a  significant  decrease  in  the 
accuracy  of  the  predicted  IMF. 

During  these  13  rotations  the  computed  current  sheet 
had  a  quasi-stationary  structure  with  two  northward  and  two 
southward  excursions  per  rotation,  corresponding  to  a  four 
sector  structure.  Occasionally  an  excursion  "missed"  the 
Earth.  Near  sunspot  minimum  the  maximum  extent  in  latitude 
of  the  current  sheet  was  about  15  degrees,  but  1.5  years 
later  the  maximum  latitude  had  increased  to  about  45 
degrees . 

Comparisons  with  the  IMF  polarity  observed  at  Earth 
give  only  a  weak  test  of  the  extent  in  latitude  of  the 
current  sheet.  The  most  definitive  answer  to  this  question 
will  come  from  observations  with  spacecraft  at  larger  helio¬ 
graphic  latitudes.  Although  a  large  part  of  the  heliosphere 
is  filled  with  magnetic  flux  from  the  solar  polar  regions, 
the  structure  of  the  IMF  observed  at  Earth  is  still  closely 
related  to  the  structure  of  low-latitude  fields  on  the 


source  surface 
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Figure  Captions 

Figure  1  Computed  magnetic  field  contours  on  a  spherical 

source  surface  concentric  with  the  Sun  at  a  radius 
of  2.35  R0  for  Carrington  Rotation  1648,  beginning 
7  November  1976.  The  solid  contour  lines  represent 
field  directed  away  from  the  Sun  with  relative 
strengths  1,  5  and  10;  the  dashed  contours  represent 
field  directed  toward  the  Sun.  The  heavier  line 
near  the  solar  equator  is  where  the  direction  of 
the  computed  field  changes  from  away  to  towards, 
and  is  assumed  to  be  the  source  of  the  heliospheric 
current  sheet.  The  +•  and  -  symbols  represent  daily 
values  of  the  polarity  of  the  interplanetary  magnetic 
field  observed  at  Earth,  adjusted  for  the  five  day 
transit  time  of  solar  wind  from  Sun  to  Earth. 

Figure  2  The  same  format  as  Figure  1,  but  for  a  later  Carrington 
Rotation  1656  beginning  13  June  1977.  Note  that 
the  extent  in  latitude  of  the  computed  heliospheric 
current  sheet  extends  to  higher  latitudes  than  in 
Figure  1. 

Figure  3a  The  heliospheric  current  sheet  computed  on  a  source 
surface  at  2.35  Rg  on  nine  successive  Carrington 
Rotations,  1641-1649,  beginning  on  30  4pril  1976  to 
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4  December  1977.  Compare  Cor  example  the  current 
sheet  shown  here  for  Carrington  Rotation  1648  with 
that  shown  in  Figure  1.  Each  succeeding  base  line 
(solar  equator)  is  displaced  by  45  degrees  heliographic 
latitude.  The  +  and  -  symbols  represent  daily  values  of 
the  IMF  polarity  observed  at  Barth  allowing  for  the  five 
day  transit  time  of  solar  wind  from  Sun  to  Earth. 
Significant  disagreements  between  the  predicted  and 
observed  IMF  polarities  are  indicated  with  a  thicker 
neutral  line.  (The  first  rotation  shown  in  Figure  3a 
is  near  sunspot  minimum.) 

Figure  3b  The  same  as  Figure  3a,  but  for  the  next  nine  Carrington 
Rotations,  1650-1658,  beginning  1  January  1977  to 
7  August  1977.  Note  that  the  extent  in  latitude  of 
the  computed  current  sheet  increases  in  the  later 
rotations. 

Figure  4a  A  synoptic  map  of  the  line  of  sight  photospheric 

magnetic  field  measurements  observed  at  the  Stanford 
Solar  Observatory  for  Carrington  Rotation  1643.  The 
dates  indicate  the  central  meridian  passage  time  for 
corresponding  longitude.  The  inverted  carets  show  the 
dates  of  magnetograph  scans  which  contribute  to  the 
chart. 


Figure  4b 


The  same  format  as  figure  4a  for  Carrington  Rotation 
1644.  Notice  the  large  active  region  that  has  appeared 
in  the  northern  hemisphere  near  longitude  120. 


Figure  4c  Synoptic  chart  for  Carrington  Rotation  1645.  The  size 
and  strength  of  the  active  region  is  greatly  reduced. 

Figure  5  Maximum  correlation  between  the  IMF  polarity  predicted 

from  the  computed  heliospheric  current  sheet  and  the 
IMF  polarity  observed  at  Earth  as  a  function  of  the 
source  surface  radius  on  which  the  current  sheet  was 
computed.  Source  surfaces  were  computed  with  an  added 
solar  polar  field  strength  of  11.5  gauss  as  computed  by 
Svalgaard  et  al.  (1978),  and  for  other  values  of  the 
added  solar  polar  field  as  shown.  , 

Figure  6  Computed  heliospheric  current  sheets  on  Carrington 
Rotation  1656  beginning  13  June  1977  for  several 
values  of  added  solar  polar  magnetic  field.  As  the 
strength  of  the  polar  field  Is  increased  the 
computed  current  sheet  approaches  the  plane  of  the 
solar  equator. 

Figure  7  Computed  current  sheets  for  Carrington  Rotation  1656 
beginning  13  June  1977  for  source  surfaces  at  several 
different  radii,  as  indicated.  As  the  radius  of  the 
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Figure 


Figure 


Figure 


source  surface  is  increased  the  computed  current  sheet 
approaches  the  solar  equator. 

8  The  maximum  cross  correlation  between  the  IMF  polarity 
predicted  from  a  computed  current  sheet  on  a  source 
surface  at  2.35  R#  with  11.5  gauss  added  polar  field 
and  the  IMF  observed  at  Barth  as  a  function  of  the 
latitude  on  the  source  surface  at  which  the  field 
polarity  was  predicted.  In  the  abscissa  zero  represents 
the  heliographic  latitude  of  the  earth. 

9  Cross  correlation  between  the  IMF  polarity  predicted 
from  the  adopted  computation  of  the  heliospheric  current 
sheet  and  the  polarity  observed  at  Barth.  The  lag  of 
the  first  peak  is  five  daysr  which  represents  the  transit 
time  from  Sun  to  Earth  of  the  solar  wind  near  the  sector 
boundaries. 

10  The  IMF  polarity  computed  at  the  source  surface  by  the 
model  is  presented  in  the  Bartels  chart  format.  Bach 
row  has  27  boxes  with  the  polarity  for  each 

day  indicated  in  a  box.  K  filled  box  indicates  toward 
polarity;  a  hatched  box  indicates  indeterminate  polarity; 
an  empty  box  indicates  away  polarity.  The  plot  is  displaced 
by  five  days  to  account  for  the  solar  wind  transit  time 
from  Sun  to  Barth.  This  format  emphasizes  the  27-day 
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recurrence  pattern  in  the  polarity  and  the  large-scale 
structure  over  many  rotations. 

Pigure  11  Same  format  as  Pigure  10,  but  for  the  IMP  polarity 
observed  at  Barth. 
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CORRELATION  OF  PREDICTED  AND  OBSERVED  IMF  POLARITY 
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